Self-organized and driven phase synchronization in coupled maps 
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We study the phase synchronization and cluster formation in coupled maps on different networks. 
We identify two different mechanisms of cluster formation; (a) Self-organized phase synchronization 
which leads to clusters with dominant intra-cluster couplings and (b) driven phase synchronization 
which leads to clusters with dominant inter-cluster couplings. In the novel driven synchronization 
the nodes of one cluster are driven by those of the others. We also discuss the dynamical origin of 
these two mechanisms for small networks with two and three nodes. 
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Recently, there is considerable interest in complex sys- 
tems described by networks or graphs with complex 
topology [1]. Most networks in the real world consist 
of dynamical elements interacting with each other. Thus 
in order to understand properties of such dynamically 
evolving networks, we study a coupled map model of dif- 
ferent networks. Coupled maps show rich phenomenology 
that arises when opposing tendencies compete; the non- 
linear dynamics of the maps which in the chaotic regime 
tends to separate the orbits of different elements, and the 
couplings that tend to synchronize them. Coupled map 
lattices with nearest neighbor or short range interactions 
show interesting spatio-temporal patterns, and intermit- 
tent behavior [2]. Globally coupled maps (GCM) where 
each node is connected with all other nodes, show in- 
teresting synchronized behavior [3]. Ref. [4] are some 
of the papers which shed light on the collective behav- 
ior and synchronization of coupled maps/oscillators with 
local and non-local connections on different networks. 

In this paper we study the mechanisms for synchroniza- 
tion behaviour of coupled maps on different networks. In 
particular, we concentrate on networks with small num- 
ber of connections, i.e. the number of connections (N c ) 
are of the order of the number of nodes (N) . Our study 
reveals two different ways for the formation of synchro- 
nized clusters, (a) Synchronized clusters can be formed 
because of intra-cluster couplings. We will refer to this as 
self-organized synchronization, (b) Synchronized clusters 
can be formed because of inter-cluster couplings. Here 
nodes of one cluster are driven by those of the others. 
We will refer to this as driven synchronization. We are 
able to identify ideal clusters of both types, as well as 
clusters of the mixed type where both ways of synchro- 
nization contribute to cluster formation. We will discuss 
several examples to illustrate both types of clusters. Dy- 
namically, our analysis indicates that the self-organized 
behaviour has its origin in the decay term arising due 



to intra-cluster couplings in the dynamics of the differ- 
ence variables while the driven behaviour has its origin 
in the cancellation of the inter-cluster couplings in the 
dynamics of the difference variables. 

Consider a network of N nodes and iV c connections. 
Let each node of the network be assigned a dynamical 
variable x l , i = 1, 2, . . . , N. The evolution of the dynam- 
ical variables is given by 



H+i 



= {l-e)f{x\) + 



where x l n is the dynamical variable of the i-th node at 
the t-th. time step, C is the adjacency matrix with ele- 
ments Cij taking values 1 or depending upon whether 
i and j are connected or not. The matrix C is symmetric 
with diagonal elements zero. The function f(x) defines 
the local nonlinear map and the function g(x) defines the 
nature of coupling between the nodes. In this paper, we 
present the results for the local dynamics given by the 
logistic map f(x) = /ix(l — x) and two types of coupling 
functions, (i) g(x) = x and (ii) g(x) — f(x). 

Synchronization of coupled dynamical systems may be 
denned in various ways [5] . Perfect synchronization cor- 
responds to the dynamical variables for different nodes 
having identical values. Phase synchronization corre- 
sponds to the dynamical variables for different nodes 
having values with some definite relation [6]. For net- 
works with N c ~ N, we find that perfect synchronization 
leads to clusters with very small number of nodes, while 
phase synchronization gives clusters with large number 
of nodes. Here, we concentrate on phase synchronized 
clusters. We define the phase synchronization as follows 
[7]. Let rij and rij denote the number of times the vari- 
ables x \ and x 3 t , t — 1,2, . . . ,T for the nodes i and j show 
local minima during the time interval T. Let fiy denote 
the number of times these local minima match with each 
other. We define the phase distance between the nodes 
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i and j as dij = 1 — 2riij/(ni + nj). Clearly, dij = formation we define two quantities j 'inter and f intra as 
when all minima of variables x % and x° match with each 
other and da = 1 when none of the minima match. We 



say that nodes i and j are phase synchronized if di 



0. 



Also, a cluster of nodes is phase synchronized if all pairs 
of nodes of that cluster are phase synchronized. 

We find examples of both self-organized and driven 
types of phase synchronized clusters in different networks 
that we have studied. For small coupling strengths, we 
observe turbulent behaviour, i.e. no clusters are formed, 
but as the coupling strength increases phase synchro- 
nized clusters are formed. The number and sizes of 
clusters as well as their type (self-organized, driven or 
mixed) depends on the coupling strength e as well as 
the type of coupling function g(x). For networks with 
number of connections of the order of N, and for linear 
coupling g(x) = x, we observe self-organized phase syn- 
chronized clusters for small coupling strengths (e ~ 0.18) 
and driven phase synchronized clusters for large coupling 
strength with a crossover and reorganization of nodes be- 
tween the two types as e is increased. This behaviour 
appears to be approximately independent of the type 
of network. On the other hand, for nonlinear coupling 
g{x) = fix), we observe a dominant driven phase syn- 
chronization. In this case, the sizes and number of clus- 
ters depends on the type of network for large e values. 
As noted earlier, in this letter we concentrate only on the 
mechanism of cluster formation and other details will be 
discussed elsewhere [8] . 

We now present the numerical results of our model. 
Starting from random initial conditions and after an ini- 
tial transient, we study the dynamics of Eq. (1), and 
determine synchronization behaviour. Fig. 1 shows sev- 
eral examples of clusters illustrating different behav- 
iors. Fig. 1(a) shows two clusters with an ideal self- 
organized phase synchronization. We note that except 
one coupling, which must be present since our network 
is connected, all other couplings arc of intra-cluster type. 
Fig. 1(b) shows the opposite behaviour of two clusters 
with an ideal driven phase synchronization. Here, all 
the couplings are of the inter-cluster type. Fig. 1(c) to 
1(e) show mixed behaviour. Fig. 1(c) shows clusters of 
different types. The largest two clusters have approx- 
imately equal number of inter-cluster and intra-cluster 
couplings (mixed type), the next two clusters have dom- 
inant intra-cluster couplings (self-organized type) while 
the remaining three clusters have dominant inter-cluster 
couplings (driven type). Also there are several isolated 
nodes. Fig. 1(d) shows clusters where driven behaviour 
dominates. Fig. 1(e) shows clusters where self-organized 
behavior dominates. Fig. 1(f) shows two clusters of ideal 
driven type with several isolated nodes. Figs. 1(c) to 1(f) 
have isolated nodes which do not belong to any clus- 
ter. These nodes evolve independently, however, some of 
them can get attached to some clusters intermittently. 

To get a quantitative picture of the two ways of cluster 



fintra 
/inter = 



Ni, 



N c 
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(2a) 
(2b) 



where N in tra is the number of intra-cluster couplings and 
N i n ter is the number of inter-cluster couplings. Couplings 
between isolated nodes are not counted in A T inter- 

Figs. 2(a) and 2(b) show both j intra and f inter as a 
function of e for g(x) — x and g(x) = f(x) respectively 
for the scale-free networks. ^From Fig. 2(a) (g(x) = x), 
we see that after an initial turbulent phase we get clus- 
ters with large values of fintra, i-e. self-organized clus- 
ters, for e > 0.12. fintra becomes almost one for e ~ 0.18 

and then starts decreasing. As e increases further f inter 
starts increasing and there is a crossover and reorganiza- 
tion of nodes to driven clusters, so that for very large e, 
f inter is close to one. The ideal driven cluster shows 
two driven clusters which are anti-phase synchronized 
with each other. On the other hand, for g(x) — f{x), 
Fig. 2(b) shows that j 'inter > i-e. driven behaviour, domi- 
nates. Fig. 2(c) and 2(d) show similar graphs for network 
with one dimensional nearest neighbor couplings. The 
behaviour is similar to that of scale-free networks except 
that for g(x) = f(x) and for large e there is almost no 
synchronization or cluster formation (Fig. 2(d)). 

It is interesting to note that the two different ways of 
cluster formation are observed even when the variables 
in the clusters are evolving chaotically. For ji = 4, we 
find that when three or more clusters are formed the 
largest Lyapunov exponent is positive. When two clus- 
ters (with or without some isolated nodes) are formed 
largest Lyapunov exponent can be both positive or neg- 
ative depending on the parameter values and the type of 
coupling. If the largest Lyapunov exponent is negative, 
the variables show periodic behavior with even period [8] . 
For mil < 4, we find different periods including odd ones 
and also two, three or more stable clusters depending 
upon the parameters and initial conditions. 

Geometrically, the organization of the network into 
couplings of both self-organized and driven types is easy 
to understand for the networks with tree structure. A 
tree can always be broken into two or more disjoint clus- 
ters with only intra-cluster couplings by breaking one or 
more connections. Clearly, this splitting is not unique. 
A tree can also be divided into two clusters by putting 
connected nodes into different clusters. This division is 
unique and leads to two clusters with only inter-cluster 
couplings. For other types of networks splitting into clus- 
ters with ideal intra-cluster or inter-cluster couplings may 
not be always possible. However, clusters with dominant 
couplings of either intra-cluster or inter-cluster type is 
still possible for N c ~ N. For larger values of N c , typi- 
cally of the order of A 2 , a clear identification of only one 
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type of behavior becomes difficult and the clusters are 
mostly of the mixed type. Note that geometrically it is 
always possible to get one big cluster spanning almost all 
the nodes of the self-organized type. 

A comment on the choice of T used to determine the 
phase synchronization. Clearly, T should be large enough 
to include several maxima and minima. On the other 
hand it should be small enough to include the behavior 
of some isolated nodes that get attached to some clusters 
intermittently with a time scale of t s . We find that r s is 
about few thousand iterates. Hence, we chose T = 100. 

We have also studied networks with large N (largest N 
was 10,000). We can clearly identify both self-organized 
and driven behavior in such large networks also. 

To understand the dynamical origin of the self- 
organized and driven phase synchronization let us first 
consider a network of two variables, x 1 and x 2 , coupled to 
each other. Synchronization between these two variables 
will be decided by the difference variable x s ~ = x 1 — x 2 . 
From Eq. (1) the dynamics of x s ~ is given by 

= (1 - e)(f{x\) f(x 2 )) 6 -( g (x\) - g(x 2 ). (3) 

Ref. [5] discusses synchronization properties of this net- 
work of two variables for g(x) = f(x). Next consider 
a network of three variables with both x 1 and x 2 cou- 
pled to x 3 and no coupling between x 1 and x 2 . Now the 
dynamics of the difference x d ~ = x 1 — x 2 is given by 

x d t ^ = (l-e){f{xl)-f{x 2 )). (4) 

It can be shown that there is a critical value of e above 
which the variables x 1 and x 2 will synchronize, i.e. x d ~ 
will tend to zero. The detailed dynamics of the above 
two simple networks and their synchronization behavior 
will be discussed elsewhere [8] . 

Comparison of Eqs. (3) and (4) clearly shows the dif- 
ferent dynamical origins of the self-organized and driven 
mechanisms. The intra-cluster coupling term which is 
responsible for the self-organized behaviour, adds a de- 
cay term to the dynamics of x s ~ (Eq. (3)). On the other 
hand, the inter-cluster coupling terms, which are respon- 
sible for the driven behaviour, cancel out and do not add 
any term to the dynamics of x d ~ (Eq. (4)). We feel that 
for larger networks also similar mechanisms as in Eqs. (3) 
and (4) are responsible for the cluster formation of the 
self-organized and driven type. 

There are several examples of self-organized and driven 
behavior in naturally occurring systems. Self-organized 
behavior is common and is easily observed. Examples 
are social, ethnic and religious groups, political groups, 
cartel of industries and countries, herds of animals and 
flocks of birds, different dynamical transitions such as 
self-organized criticality etc. The driven behavior is prob- 
ably not so common [10]. An interesting example is the 
behavior of fans during a match between traditional ri- 
vals. Before the match the fans may act as individuals 



(turbulent behavior) . During the match, when the match 
reaches a feverish pitch, i.e. the strength of the inter- 
action increases, the fans are likely to form two driven 
phase synchronized groups. The response of each group 
depends on that of the other and is normally anti-phase 
synchronized with the other. Another example is the for- 
mation of opposite ethnic groups as in Bosnia. 

In this letter we have presented results for the case 
where the local dynamics is governed by logistic map 
and couplings g{x) — x and f(x) for the scale free net- 
works and one-d lattice with nearest neighbor couplings. 
We have also studied several other maps (circle map, 
tent map etc.), other types of couplings and different 
networks (two-d lattice with nearest neighbor couplings, 
small world networks, Caley tree, random networks, bi- 
partite networks) . We find similar behaviour and are able 
to identify self-organized and driven behaviour. 

To conclude we have investigated the mechanism of 
cluster formation in coupled maps on different networks. 
We are able to identify two distinct ways of cluster for- 
mation, namely self-organized and driven phase synchro- 
nization. Self-organized synchronization is characterized 
by dominant intra-cluster couplings while driven behav- 
ior is characterized by inter-cluster couplings. The two 
ways of cluster formation are clearly seen for networks 
with small number of couplings (N c ~ N) but are diffi- 
cult to identify as the number of couplings increases and 
becomes of the order of A 2 . Dynamically, the examples 
of small networks show that the self-organized behaviour 
occurs because of the intra-cluster couplings introducing 
a decay term in the difference variables while the driven 
behaviour occurs because of the inter-cluster couplings 
cancelling out. 

One of us (REA) thanks Professor H. Kanz for useful 
discussions and Max-Plank Institute for the Physics of 
Complex Systems, Dresden for hospitality. 
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FIG. 1. The figure shows several examples illustrating the self-organized and driven phase synchronization. The examples are 
chosen to demonstrate the two different ways of obtaining synchronized clusters and the variety of clusters that are formed. All 
the figures show node verses node diagram for N = Nc = 50. After an initial transient (about 2000 iterates) phase synchronized 
clusters are studied for T = 100. The logistic map parameter fi = 4. The solid circles show that the two corresponding nodes 
are coupled and the open circles show that the corresponding nodes are phase synchronized. In each case the node numbers are 
reorganized so that nodes belonging to the same cluster are numbered consecutively, (a) Figure show an ideal self-organized 
phase synchronization for scale free network for g(x) = x ande = 0.15. (b) An ideal driven phase synchronization for scale free 
network for g(x) = x and e = 0.85. (c) Mixed behavior for icale free network for g(x) = f(x) and e = 0.61. (d) A dominant 
driven behavior for scale free network for g(x) — f(x) and e = 0.87. (e) A dominant self-organized behavior for 1-d lattice with 
nearest neighbor couplings for g(x) — x and e = 0.14. (f) An ideal driven behavior with several isolated nodes for 1-d lattice 
with nearest neighbor couplings for g(x) — f(x) and e = 0.15. The scale free networks were generated starting with No = 1 
nodes and adding one node with m = 1 couplings at each stage of the growth of the lattice with probability of connecting to a 
node being proportional to the degree of the node (see Ref. [9] for details). 
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FIG. 2. The fraction of intra-cluster and inter-cluster couplings, fi„ter (solid line) and fintra (dashed line) are shown as a 
function of the coupling strength e. Figures (a) and (b) are for the scale free network for g(x) = x and g(x) = f(x) respectively. 
Figures (c) and (d) are for the 1-d network with nearest neighbor couplings for g(x) = x and g(x) = f(x) respectively. The 
figures are obtained by averaging over 20 realizations of a network and 50 random initial conditions for each realization. 
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